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Abstract 

We propose a cosmological model which exhibits the phenomenon of self-accele- 
ration: the Universe is attracted to the phase of accelerated expansion at late times 
Qh' even in the absence of the cosmological constant. The self-acceleration is inevitable in 

Q : 

the sense that it cannot be neutralized by any negative explicit cosmological constant. 
The model is formulated in the framework of brane-world theories with a warped extra 
dimension. The key ingredient of the model is the brane-bulk energy transfer which 
is carried by bulk vector fields with a a-model-like boundary condition on the brane. 
We explicitly find the 5-dimensional metric corresponding to the late-time de Sitter 
expansion on the brane; this metric describes an AdSs black hole with growing mass. 
^^ ■ The present value of the Hubble parameter implies the scale of new physics of order 
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The Universe expansion history predicted by the model is distinguishable from the 
standard ACDM cosmology. 



1 TeV, where the proposed model has to be replaced by putative UV-completion. The 
mechanism leading to the self-acceleration has AdS/CFT interpretation as occurring 
due to specific dynamics of conformal matter interacting with external "electric" fields. 



1 Introduction and summary 

One of the open questions of the contemporary physics is the nature of dark energy. The 
latter term is conventionally used for the unknown substance which is believed to be respon- 
sible for the observed acceleration of the cosmological expansion. In the minimal scenario, 
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the dark energy is identified with the cosmological constant (CC). However, the energy scale 
associated with the value of the CC required to explain the observed dark energy density [1], 

p DE = (0.43 • fCT 10 GeV) 4 , (f ) 

is much smaller than any other fundamental energy scale existing in particle physics. Thus, 
one faces the problem of explaining this strong hierarchy [2]. All particle condensates (i.e., 
those associated with spontaneous symmetry breaking) contribute to CC but have to almost 
cancel out. Moreover, the CC receives contributions from quantum fluctuations of all fields 
so that a natural value for it would be 

^natural) „ ^ (2) 

where hjjv is the ultra-violet cutoff of the theory. The validity of the Standard Model of 
particle physics up to energies of order 1 TeV implies that Ajjv <^ 1 TeV. The latter value of 
Kuv when substituted into (2) gives the natural value of the CC much larger than the one 
consistent with observations. 

Of course, Eq. (2) is not rigorous: quantum corrections to CC are divergent and can be 
formally put to any value by renormalization. There exist several proposals to cancel the CC 
(for reviews see, e.g., Refs. [2, 3, 4]). They can be split in two groups. One group involves 
proposals based on relaxation mechanisms [5, 6, 7, 8, 9]. These proposals suffer from the 
"empty Universe" problem (see, e.g., [4]): the relaxation mechanisms drive the Universe to 
the state practically devoid of matter which is inconsistent with the Hot Big Bang. It is not 
clear for the moment whether this problem can be consistently resolved without invoking 
anthropic selection [10]. Another group involves symmetry [11] or quantum gravity [12, 13] 
arguments. However, these proposals adjust the CC somewhat too well: they predict that 
the CC is exactly zero. The problem is now reversed: obtaining the Universe undergoing the 
late-time de Sitter expansion requires fine-tuning. 

In this context one would like to construct a setup which excludes the Minkowski space- 
time from the possible asymptotics of the history of the Universe irrespectively of the value 
of the CC. Namely, we envisage the following picture. The Universe is driven to the de 
Sitter phase at late times not by the CC but by some dynamical mechanism. The Hubble 
parameter of the asymptotic de Sitter space-time may depend on the value of the CC but is 
always bounded from below by a tiny positive value. This reduces the CC problem to its old 
form: the need to explain the cancellation of the CC (admittedly, this is the most difficult 
task) without the risk of being over-precise. An additional requirement one would like to 



impose on the dynamical mechanism leading to the late-time acceleration of the Universe is 
that it should not invoke the small energy scale entering Eq. (1) as an input and should not 
require any particular initial conditions for the evolution of the Universe. It is natural to 
refer to the scenario described above as "self-accelerated cosmology". 

The model that came closest to the realization of the self-accelerated scenario is the DGP 
model [14]. It is formulated in flat (4+l)-dimensional space-time with our Universe being 
represented by a 3-dimensional brane; the key ingredient of the model is the induced gravity 
term on the brane. Still, the self-accelerated branch of this model suffers from the presence 
of ghost fields [15, 16], i.e. excitations with negative kinetic term which signal inconsistency 
of the theory. 

In this paper we propose another model which realizes the idea of self-accelerated cos- 
mology. The model is formulated in the context of brane-world scenario with warped extra 
dimensions [17]. We do not consider the induced gravity term on the brane. Instead, we 
introduce an energy exchange between the brane and the bulk carried by a triplet of U(l) 
bulk vector fields [18]. Then, the self-acceleration mechanism can be summarized as follows. 
As the brane expands, it gives rise to an energy flow from the brane into the bulk. This 
energy flow, in its turn, results in the formation of a bulk black hole with growing mass. 
The latter back-reacts on the brane and maintains the brane Hubble parameter constant at 
late times. 3 Similar ideas were expressed previously in Ref. [20]; our model provides explicit 
realization of these ideas in the framework of field theory. The brane-bulk energy exchange 
is kept at the level sufficient to provide the self-accelerated cosmology by a specific Dirichlet 
boundary condition for the vector fields which forces their squared values on the brane to be 
equal to a fixed parameter M 2 , see Eq. (5). Note that existence of a continuous energy flux 
from the brane implies violation of the null energy condition by the brane energy-momentum 
tensor. This property is dangerous as it can lead to instabilities [21]. Below we briefly dis- 
cuss this problem and a possible way to resolve it. A more detailed treatment of this issue 
will be presented elsewhere [22]. In the present paper we concentrate on the cosmological 
applications of our model. 

When the explicit CC on the brane vanishes, the self-accelerated expansion regime at 
late times corresponds to the value of the would-be CC given by the seesaw formula, 
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(M 2 /M P/ ) 4 . (3) 



3 The increase of the brane Hubble parameter by a bulk black hole is a known effect: a static black hole 
acts as "dark radiation" term in the Friedman equation [19]. For a constant mass black hole this effect dies 
out at late times while in our case it persists due to the growth of the black hole mass. 



The present dark energy density (1) corresponds to M of order 1 TeV. It is important to 
stress that self-acceleration mechanism we propose is not equivalent to an implicit introduc- 
tion of a small positive CC. We check this in two ways. First, we find that in the absence of 
an explicit CC and matter on the brane the system admits, apart from the self-accelerated 
solution, also a static solution with Minkowski metric on the brane. However, in contrast 
to the self-accelerated solution, the Minkowski solution is unstable and is not an attractor 
of the cosmological evolution. This solves the problem of naturalness of the initial condi- 
tions leading to late-time acceleration: in our model the Universe always approaches the de 
Sitter space at late times. Second, we show that the self-acceleration mechanism cannot be 
compensated by an explicit introduction of a negative CC on the brane. Namely, depending 
on the value of the negative CC, the Universe either undergoes the accelerated expansion 
with essentially the same Hubble parameter corresponding to (3), or collapses. The lat- 
ter behavior, the collapse for large and negative explicit CC, is similar to the collapse of 
the Universe in the case of negative CC in the standard Friedman cosmology. Presumably, 
a (yet unknown) mechanism which cancels the large positive CC, excludes the collapsing 
cosmological solutions as well, cf. [12, 13]. 

The self-acceleration mechanism realized in our model has a transparent interpretation 
in term of the AdS/CFT correspondence [23]. In this language the acceleration is driven by 
the conformal matter. The cooling of this matter due to the expansion of the Universe is 
compensated by an energy inflow due to work done by external "electric" fields on the CFT. 
The latter do not decay rapidly in the expanding Universe in our model. Technically, this 
is achieved by the specific a-model-like conditions on the vector fields, Eq. (5). The CFT 
interpretation of the self-acceleration mechanism suggests phenomenological generalizations 
of the model. The cosmological expansion equation in this class of models is distinct from 
the standard ACDM case, thus providing falsifiable predictions. A careful analysis of this 
and other phenomenological consequences of our model and its generalizations is left for 
future work. 

The paper is organized as follows. In Sec. 2 we introduce the setup and identify its 
features that are relevant for cosmology. Section 3 contains qualitative analysis of the cos- 
mological evolution in the model. A more rigorous analysis is performed in Sec. 4 where 
an explicit solution describing a self-accelerated universe is given. In Sec. 5 we describe the 
interpretation of this cosmological solution in the spirit of the AdS/CFT correspondence. In 
Sec. 6 we discuss open issues related to our work. Appendix A presents an explicit example 
of a brane Lagrangian giving rise to the self-accelerated cosmology. Appendices B and C 



contain technical details. 

2 Setup 

We consider a brane- world model based on the setup analogous to that of Ref. [17]. All the 
fields describing the observable matter are supposed to be localized on a 3-brane. The action 
of the model is the sum of bulk and brane parts, 

The bulk sector consists of gravity described by the standard Einstein-Hilbert action with 
negative cosmological constant A 5 and of a triplet of Abelian bulk vector fields A a M) a = 
1, 2, 3. Thus, for the bulk action we have, 4 
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where 

Fmn = OmA n — OnA m . 

The brane action is taken in the following form, 

Sbrane = / d^Xy^j ( - O + C S M + ^Hs[^ •••])• 



(4) 



Following [17] we impose Z 2 reflection symmetry across the brane. In the above brane action 
a is the tension of the brane, Csm represents the Standard Model Lagrangian (possibly, 
including dark matter sector), and Chs is a boundary Lagrangian for the vector fields. The 
dots represent possible dependence of Chs on some additional fields localized on the brane. 
We will refer to the sector described by Chs as the "hidden sector": we assume that there 
are no direct interactions of the vector fields A a as well as of the other fields entering into 
Chs (except the metric g^ u ) with the Standard Model fields. 
For the brane tension we write 
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v /-3A 5 /AnG 5 + p A . 



4 Our convention for the metric signature is (+, —,—,...); we use upper-case Latin letters M, N, . . . for 
5-dimensional indices, the Greek letters /i, v, . . . for 4-dimcnsional indices in the tangent space to the brane, 
and lower-case Latin letters i,j,... for the 3-dimcnsional spatial indices. 



As we will see shortly p\ plays the role of the 4-dimensional CC. We assume that there 
is some mechanism ensuring that p\ is much smaller than its natural value (2). Without 
knowledge of the details of this mechanism we cannot be more concrete and say what value 
of Pa h predicts. A reasonable possibility is that p\ is set equal to zero but other options 
cannot be excluded. So we will keep p\ explicitly in the formulas in order to be able to 
control its effect on the cosmological evolution. 

The precise form of the hidden sector Lagrangian is not important for our purposes. For 
the study of cosmology it is sufficient to require the following properties: 

(1) Variation of Chs gives the following boundary condition for the vector fields on the 
brane: 

g^A a A b u = -M 2 5 ab , (5) 

brane 

where g^ v is the induced metric on the brane, and M is some parameter of dimension of 
mass. Though the condition (5) has a generally covariant form it can be satisfied only if the 
vector fields are non-zero so that the local Lorentz symmetry is spontaneously broken. Note 
also that Eq. (5) explicitly breaks the U(l) gauge symmetries of the bulk action. Having in 
mind an analogy with cr-models of particle physics one expects that the use of the condition 
(5) is justified at energies below M which are relevant for late-time cosmology. On the other 
hand, at energies above M some UV-completion of the theory is required. 5 

(2) The energy density of the hidden sector matter dilutes, pus — ¥ 0, as the brane 
undergoes the cosmological expansion. In other words, we assume that the hidden sector as 
well as the Standard Model sector does not contain any CC. In order not to contradict the 
observational data pns should dilute at the same rate or faster than the energy density of 
dust. 

An explicit example of the Lagrangian Chs satisfying the requirements (1), (2) is pre- 
sented in appendix A. 

We make the following assumptions about the order-of-magnitude values of the parame- 
ters of our setup. The bulk physics is governed by the energy scale of order or slightly below 
the Planck mass. Thus, we have 

G 5 ~ M~f , e 2 - M" 1 , A 5 - M 5 Pl . (6) 

On the other hand, the scale M entering into Eq. (5) is assumed to be much smaller than 
the Planck mass. Below we will take M ~ 1 TeV in order to fit the observed acceleration 



5 In fact, the scale M is an upper bound on the UV-cutoff of the theory. The actual value of the cutoff can 
be lower in a given model of the brane Lagrangian. This issue goes beyond the scope of the present article. 
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rate of the Universe. 

3 Cosmology: qualitative discussion 

3.1 Cosmological Ansatz 

We now turn to the discussion of cosmology in the model described above. The general 
Ansatz for the metric and vector fields, invariant under 3-dimensional translation and rota- 
tions, has the form 

ds 2 = F{tX){dt 2 - d( 2 ) - r 2 (t,()(dxf; (7) 

A* = MA(t, C) 8° , A a t = A a c = 0. (8) 

Equations (8) imply that the vector fields A a M form an orthogonal spacelike triad; such 
a configuration is invariant under 3-dimensional rotations supplemented by simultaneous 
rotations in the internal space. Embedding of the brane into the space with metric (7) is 
described by a set of functions x M (y fJ- ). Choosing x % = y l one is left with two functions, 
£b(t) and (b{t), describing the trajectory of the brane on the (t, C) plane. Here r is a time 
parameter on the brane, see Fig. 1. It is convenient to choose r to be the cosmological time. 
Then the induced metric on the brane takes the Robert son- Walker form 

ds 2 = dT 2 -a 2 (r)(dx 1 ) 2 . (9) 

Here a(r) = r(is(7"), Cb(t)), and the functions ts, (b obey the condition 

F(t B ,(B){i 2 B -( 2 B ) = l, 

where dot denotes derivative with respect to r. The motion of the brane is affected by the 
bulk metric via the junction conditions. Imposing Z 2 reflection symmetry across the brane 
one obtains, 

47rG 5 T ( ; r)T = -3 (^( B + ^t B ) , (10a) 
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2d t r d t F\ ■ (28 c r d c F\ . 1 d r-, 

— + -2-F-) CB + {— + -W) tB + 7fTb^ Cb) 



sj. (iob) 



Finally, the boundary conditions (5) relate the boundary value of the function A(t, () to the 
scale factor of the brane metric, 

A(t B (r)XB(r))=a(r). (11) 



tk 




Figure 1: The coordinate frame used in the text. Thick curve shows the trajectory of the 
brane. The space-time considered in the text consists of two copies of white region (above 
the thick curve) glued across the brane trajectory and related by Z2 symmetry Gray region 
is cut off. 



This equation implies that the cosmological expansion of the brane acts as a source for the 
bulk vector fields. As we will see, this lies at the origin of the self-acceleration mechanism 
in our model. 

The junction conditions (10) are simplified if the Hubble parameter H = a/ a on the 
brane is small compared to the curvature of the bulk. In this case one expects the distortion 
of the bulk due to the vector fields emitted by the brane to be small near the brane. Then 
the metric in the vicinity of the brane is approximately that of AdSs: 



l + 20(t,Q 



mo 



k( 



(12) 



where k = y^—AnG^A 5 /3 is the inverse AdS radius, and (j), ip "C 1. Expanding Eqs. (10) to 
the linear order in 0, ip and to quadratic order in H/k one finds the equations governing the 



cosmo logical expansion of the brane, 

H 2 = ^(PA + P) + 2A; 2 (0 + W)L ne , (13) 

H = -AnG N (p + p) + k 2 (Cd ( <f> - (d^)\ brane . (14) 

Here GV = G5A; is the four-dimensional Newton constant, and p, p are the energy density 
and pressure of the brane matter (including both the Standard Model and hidden sector 
fields). Without the terms on the r.h.s. containing the functions 0, ip, equations (13), (14) 
would describe the standard cosmology with the cosmological constant p\. The above terms 
are responsible for the modification of the brane cosmology in our model; they are to be 
found from the bulk dynamics. 

Let us comment on the validity of the approximation (12). We will check below that it 
is indeed valid in the vicinity of the brane if 

H/k < 1 . 

The latter inequality is readily satisfied in the regime of interest: eventually we want H to 
be comparable to the present Hubble parameter of the Universe, while k is assumed to be 
of order Mp\. On the other hand the form (12) is inapplicable far from the brane: we will 
see shortly that the effect of the vector fields on the space-time geometry increases deep into 
the bulk and leads to strong deviation of the metric from that of AdSs. 

3.2 Vector fields in external metric 

Let us analyze the behavior of the vector fields in the curved metric (7). Using the cosmo- 
logical Ansatz (8) one derives their equation of motion 

d t (rd t A) - d c (rd c A) = (15) 

and the expressions for the components of their energy-momentum tensor: 

Tu = T a = —((d t Af + (d c Af), (16a) 

3M 2 
T K = ^d t Ad,A , (16b) 

M 2 
Tn = S l3l ^((d t A) 2 - (d c A) 2 ) . (16c) 

Note that non- vanishing value of T t ( implies non- vanishing energy flux along the ^-direction 
carried by the vector fields. 

9 



It is instructive to neglect for the moment the backreation of the vector fields on the 
metric and consider their evolution in the external AdSs geometry. Accordingly, one sets 

1 1 

(17) 



Let us also forget about the brane for a while and consider the Dirichlet boundary condition 
for the vector field at the AdS boundary, 

A(t,()^AM(t), C^0+. (18) 

The relevance of this problem for the original cosmological setup will become clear shortly. 
Performing the Fourier decomposition, 



iujt 



A(t,() = J^A(tu,()e- 
and solving Eq. (15) one finds 

i(c,o = ycCi/ 1 (1) K)i (0) M- (19) 

This expression suggests two important observations. 

First, the vector fields are almost constant at distances ( <C 1/uj from the AdS boundary, 

^,0-^,0) «1, for. C «l. 
A(u,0) 

The latter condition is satisfied by the position (b of the brane in our cosmological setup. 
Indeed, the characteristic frequency entering the problem can be estimated as 

^ a( \" A{ \ bXb) . (20) 

A(t B ,( B ) dt B v ; 

Using the relation (11) between the boundary value of the vector field and the scale factor 
on the brane one obtains, 

Cb da da H 

u(b ~ — -7— = Cb-t- = -r < 1 , (21) 

a dtB dr k 

where in the last equality we substituted a = 1/(/cCb)- Therefore one can simplify the set of 
equations determining the cosmological solution by replacing Eq. (11) with 

A(t B (r),0)=a(r) . (22) 

10 



Instead of being formulated on the moving brane, the boundary condition for the vector field 

is now imposed on the fixed timelike surface — the AdS boundary. 

Second, the formula (19) can be used to estimate the backreaction of the vector fields 

on the metric. For example, let us consider the trace of energy-momentum tensor of these 

fields. From Eqs. (16) one obtains that the contribution of a given harmonic at uj( ^> 1 has 

the form 

M _ toM*ku,\AW(u)\* 3 

Im ~ 64^ {K) ■ 

We see that Tjjjf grows into the bulk. Thus, at large distance from the AdS boundary 
the backreaction of the vector fields on the space-time geometry gets large and should be 
properly taken into account. We proceed to this task. 

3.3 Account for backreaction 

The above analysis suggests the following physical picture. The cosmological expansion of 
the brane generates the bulk vector fields. The T t ^ component of the energy-momentum 
tensor of these fields is non-zero, so they carry a flux of energy out of the brane. Due to 
the warping of the bulk metric the energy density of the flux grows into the bulk leading to 
formation of a black hole 6 . It is known [19] that the presence of a constant mass black hole 
in the bulk affects the cosmology on the brane via the so called "dark radiation" term which 
dilutes as 1/a 4 at late times. In our case the mass of the black hole is not constant: it grows 
with time due to the continuous inflow of energy from the brane. This leads to nontrivial 
consequences at late times. 

To determine the growth rate of the black hole mass rigorously, one needs to solve the 
system consisting of equations of motion for the vector fields and the bulk Einstein equations. 
This will be done in the next section for the case of the self-accelerated cosmological Ansatz. 
For the moment we adopt a simplified strategy which allows to catch the qualitative features 
of the resulting cosmology. Namely, we will solve for the vector fields in the background of a 
static black hole. The energy flux given by this solution will be then inserted into the energy 
conservation equation in order to calculate the dependence of the black hole mass on time. 
Finally, we will use the metric of the black hole with time dependent mass to evaluate the 
bulk terms in the effective Friedman equation (13) on the brane. 



6 Anothcr possibility could be formation of a naked singularity. According to the "cosmic censorship" 
hypothesis we assume that this is not the case. This assumption will be confirmed by explicit solution of the 
Einstein equations in Sec. 4. 
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A bulk black hole 7 with mass M. is described by the AdS-Schwarzschild metric, 



ds 2 = F(r)dt 2 - 



2 dr _ 2i 2 

¥fm r l ' 



where 



F(r) 



r - 



and the Schwarzschild radius r^ is related to the black hole mass Ai as follows, 



Th 



8G 5 k 2 M^ 1/4 



371 



(23) 



(24) 



(25) 



This metric can be cast into the form (7) with ( determined by the differential equation 



Explicitly, 



c 



Akri 



In 



<K__ 

dr 


1 
" ~kF ' 






r -r h 
_r + r h _ 


1 
2kr h 


r 
arctg 

r 'h 


7[ 

' 2 



In this background the vector equation (15) takes the form 



_ d 2 A _^L drA+d 2 A={} 

<■ d( 



K' 



(26) 



(27) 



To understand the structure of the solutions to this equation let us solve it with the Dirichlet 
condition on the AdS boundary, 



One looks for the solution in the form A 
frequencies, 

uo <^ kr h . 



A(t, C = 0) = 4°) e -^ . (28) 

A^(C) e ~ lut . Let us consider the case of small 



(29) 



Then, the solution can be constructed in the following way: one solves Eq. (27) in two 
regions, ( < 1/w and ( ^> 1/kr^, and matches the solutions in the intersection of these 
regions, 1/krh «C C ^ 1/uj. At small values of ( the function A w (£) varies slowly with ( and 
one can approximate its time derivatives by their values on the AdS boundary, 



dfA^ 



-u 2 A^ . 



(30) 



7 Strictly speaking, the configuration we are considering should be called black brane as it is extended 
along 3 spatial direction. However, wc prefer the more conventional term "black hole" . 
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Equation (27) yields 



d c A„ 



4 Ax 



In 



4 



h 



+ c u 



(31) 



where C w is the integration constant. The range of validity of this formula is restricted by 
the condition (30) which translates into 

1 1 



c« 



In the region ( ^> 1/kr^ the expression (31) takes the form, 



(32) 



drA 



<^V> 



Akr h 



7[ 



31n2 -- + 0. 



U) 



! 4 0) c • 



(33) 



On the other hand, at ( ^> 1/krh the radial function r(() gets stabilized at the value of 
the Schwarzschild radius, and the second term in Eq. (27) vanishes. So, in this region the 
solutions of Eq. (27) are free waves. We take the outgoing solution 

4,(0 = A^ . 

At ( <C \juj this formula gives 

d c A w = iuAW - c/A«C . 
Comparing this expression with Eq. (33), one obtains 

A {0) = A (1) , 



Cu = -i h 3 In 2 + - w -i 

uj 2 to 



(34a) 
(34b) 



Note, that though C w is large, \CJ\ ^> 1, the second inequality in (32) is still compatible 
with ( ^> 1/krh. Substitution of (34b) into Eq. (31) yields at ( <C 1/^, 



drA 



Q^-ui 



4kr 



In 



-1 



+ 



lUJTh 



A (0) 



(35) 



Due to the condition (29) the first term on the r.h.s. of Eq. (35) is small in comparison with 
the second one if the logarithm is not too large. We will see below that this condition is 
satisfied by the solution we are interested in. 

It is straightforward to generalize the formula (35) to the case of an arbitrary dependence 
of the Dirichlet boundary condition on time, 

A(t,( = 0) = A^\t). 
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Omitting the first term on the r.h.s in Eq. (35) one obtains 

d c A = -—d t A {0) . (36) 

This expression is valid when the characteristic frequency uj ~ dfA^/A^ satisfies the in- 
equality (29). 

The next step is to determine the change of the metric due to the inflow of energy from 
the brane into the black hole. The increase of the black hole mass is given by the formula 

2tt 2 , 

d t M = -pyfiTf . 

It is convenient to evaluate the r.h.s. of this equation on the brane. Using Eqs. (16b), (22) 

and (36) one obtains 

6tt 2 M 2 2 

d * M = -fc3^ r ^°) • 

Substitution of the relation (25) yields the following time-dependence of the Schwarzschild 
radius, 

rh= {^J dt{dta)2 ) ' (37) 

where we introduced 

A = \ . 38 

Note that with our assumptions about the hierarchy among the model parameters discussed 
in Sec. 2, we have A ~ (M/M Pl ) 2 < 1. 

One proceeds by extracting the functions and ip defined in Eq. (12) from the metric 
(23) and by substituting them into the effective Friedman equation (13). One obtains, 



%nG N . . k 2 ri 



where we made use of the relation a = 1/(/c£b). Finally, Eq. (37) gives, 

H 2 = ^^( PA + p)+6, (39a) 



5=-(^ fdraa 2 ) ' . (39b) 



where 

a 4 V2A; 

In the last formula we changed the integration variable from the conformal time t to the 
proper time r on the brane using dr = adt. 
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Let us discuss our result. Equations (39) describe the cosmological evolution in terms 
of brane variables alone. The contribution of the bulk is encoded in the term (39b) which 
is nonlocal in time. This feature of our setup makes it different from other models of dark 
energy. In particular, the expansion history of the Universe at the epoch of transition from 
matter to dark energy domination is expected to be distinguishable from the standard ACDM 
case. 

The non-locality of Eq. (39a) somewhat complicates its exact solution. However, the 
qualitative effect of (39b) is easy to work out. Let us first consider the case when S is small 
and the cosmological evolution is dominated by matter with equation of state p = wp, w < 1. 
Then one has, 



p oc a~^ l+w ^ , 



and 



2 



a oc r 3 < 1+ ™) 



Substituting this expression into Eq. (39b) one notices that the integral is saturated at the 
upper limit. This implies that the choice of the integration constant is irrelevant at the late 
stages of the cosmological evolution. One obtains, 

2A \ 4/3 k 2 ( \ 2 k N 2/3 



6 \3{l-w 2 )k) r 4 /3 \(l-w) 2 Hj H ' 
We see that the 8 term dilutes slowlier than the contribution of the ordinary matter, 

5 oc a- 2 ^ . 

Thus eventually it starts to dominate the cosmological expansion. To figure out what hap- 
pens in this case we consider the opposite limit of small ordinary matter density, 8w 3 N p <C H 2 . 
The solution to Eq. (39a) is approximately given by exponential expansion, a oc e Hr and the 
Friedman equation takes the form, 

^ = ^ A + ,)+g) 2/ V (40) 

One concludes that as the matter on the brane dilutes the cosmological expansion is attracted 
to the self-accelerated fixed point. 

The asymptotic Hubble parameter depends on the explicit cosmological constant p\, see 
Fig. 2. For the case p\ = it takes the value 

H c = \ 2 k/A . 
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Pmin PA 

Figure 2: Dependence of the asymptotic Hubble parameter on the explicit cosmo logical 
constant p\. Only the solutions from the branch BC are stable and can be reached in the 
course of the cosmological evolution. 

Recalling the assumptions (6) about the parameters of the model one observes that H c 
corresponds to the effective energy density given by the seesaw formula (3). The present 
value of the Hubble parameter is obtained if M ~ 1 TeV. 

We stress that the self-acceleration mechanism described above is not equivalent to an 
implicit introduction of the CC. Indeed, for p\ = p = equations (39) admit, apart from 
the exponentially expanding solution, a static solution with H = 0. However, inspection 
of Eqs. (39) shows that this solution is unstable: addition of a small amount of matter 
drives it to the self-accelerating regime. Moreover, one cannot neutralize self-acceleration by 
introducing a negative p\. Let us discuss this issue in more detail. If the matter density p 
is zero and the cosmological constant is negative but not too large, 



> Pa > Pn 



Hi 



18ttG 



N 



equations (39) have two branches of solutions, see Fig. 2. However, one can show that 
the solutions of the branch AB are unstable: they are destroyed by the introduction of 
matter. Thus, only the solutions of the branch BC can be reached in the course of the 
cosmological evolution. The Hubble parameter on this branch is bounded from below by 
Hmin = (2/3) 3/,2 if c . If pa < p m in there are no solutions to Eqs. (39) unless there is some 
matter on the brane, p ^ 0. We solved Eqs. (39) numerically in the case of nonrelativistic 
matter and p\ < p m in and found that the Universe eventually collapses. Thus, the fate of the 
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Universe is restricted to two possibilities, depending on the value of p\: the Universe either 
undergoes a de Sitter phase with the Hubble parameter greater than H min or collapses. In 
other words, if the Universe is not collapsing it is inevitably self-accelerating. 

Let us check various approximations we made in deriving Eqs. (39). Using Eq. (37) we 
observe that the condition (29) is satisfied if 

H < V\k . (41) 

For our choice of parameters yXk ~ M and the requirement (41) is indeed satisfied for most 
stages of the cosmological evolution 8 . A problem arises with the quasistationary approxi- 
mation for the bulk metric. Indeed, from Eq. (37) one concludes that dtThjrh ~ dta/a and 
therefore the metric evolves at the same rate as the vector fields. A priori, this means that 
the formulas of this section should be taken with a grain of salt and are valid only qualita- 
tively. However, we expect the final result, Eqs. (39) to be correct. This is supported by 
the explicit solution of the time-dependent Einstein equations for the case of self-accelerated 
cosmology in Sec. 4 and the interpretation in terms of AdS/CFT correspondence in Sec. 5. 

4 Self-accelerated solution 

In this section we explicitly solve the system of bulk equations for the metric (7) with Ansatz 
(8) for the vector fields and show that the model possesses a self-accelerated solution. It is 
convenient to write the equations of motion in the light-cone coordinates 

u = k(t-C), v = k(t + 0- (42) 

A straightforward calculation yields the Einstein equations, 

d%r d v Fd v r A 2 



^(<W , (43a) 

+ ^— ~ = ~Ad u Af , (43b) 



v Ft t 

dlr d u Fd u r = A 

1" r~tl" 1" 



-d u d v In F - ^ + ^p^ = - 2 8 u Ad v A , (43c) 

d u d v r 2d u rd v r 

H — — = -F . (43d) 



The equation for the vector fields takes the form 



d u (rd v A)+d v (rd u A) = 0. (44) 



3 The inequality (41) holds when matter energy density p is smaller than (10 9 GcV) 
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One notices that the system (43), (44) is invariant under dilatations, 

A(u, v) \— > 7^(714, 7i>) , r(u, i>) i—*- jr(ju, jv) , F(u,v) 1— > 7 2 F(7«, 71;) , (45) 

where 7 is an arbitrary constant. In general, this symmetry is broken by embedding of the 

brane and by the boundary condition (11). 

An exception from this rule, which is of primary interest to us, is the case of inflating 

brane. As discussed in Sec. 3.1 the metric in the vicinity of the brane is close to the metric of 

AdSs, see Eqs. (12). The inflating brane is embedded into the AdS space along the straight 

line 

H 

C = -jt , t < , 

where we used the relation H <C &. Obviously, this embedding respects the symmetry (45). 
Due to the form of the scale factor, 

= -^, (46) 

the boundary condition (11) is also invariant under the transformations (45). 
The general Ansatz invariant under the symmetry (45) has the form 



1 (V 

A = —a 



(I), r=-I/j(l), F=l/n. (47) 

Its substitution into Eqs. (43), (44) reduces these equations to a system of ordinary differen- 
tial equations for the functions a, /3, f depending on a single variable x = v/u. One obtains, 



/3" f'(3' fa'V /A , 

-j + tH?) ' (48a) 

^_!| + ^ + 4^ + «y, (48b) 

(3 xj3 f(3 xf \(3 x(3J 

j) + J + Jf-^- A {j) =-1^) "-• (48C) 



1 


fa'V 


aa' 


1 


kJ) 


x(3 2 


/ 


1 





(3 xp \ (3 J x 

a" 13' a' 5a' (3' 
a (3a 2xa 2x/3 

The first four equations correspond to the Einstein equations, while the last one is the 
equation of motion of the vector fields; only three equations out of five are independent. 
One looks for a solution to these equations in the interval < x < 1. The end-points 
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x = and x = 1 of this interval map to the horizon ( = — t and the AdS boundary ( = 0, 
respectively. Note that we again used the inequality H <C k and identified the brane with 
the AdS boundary 9 We require that the solution should be regular at the horizon: the 
functions a, j3, f must be continuous at x — 0. The value of the function a at the AdS 
boundary is fixed by Eqs. (22), (46), (47): 



a 



a 



Vxk 



X 



(49) 



The behavior of the functions ft, f near the boundary is determined from the requirement 
that the metric in this region approximates that of AdSs. Comparing Eqs. (17) with (47) 
one obtains, 

0->t^-, f^TT^-2' *->!■ ( 5 °) 

1 — x (1 — X) 1 

Taking into account Eq. (41) we consider the case when the parameter a is large, «o ^ 1- 
The details of the solution are presented in appendix B. Here we list the results. The 
function a(x) is constant up to small corrections, a(x) m ocq. The function ft(x) is given 
implicitly by the following expression, 



2[3 h 



In 






1 



(3 n 

arctg — - - 

Ph * 



1 — X , 



where 



Ph = (ag/2) 



1/3 



Finally, the function f(x) has the form 



I = P 2 - 



Pt 

P 2 



at (1 — x) <?C 1 and 



/ oc Plx^ 



2/3 



(51) 



(52) 



(53) 



(54) 



at (1 — x) ^> 1/ 'Ph- The last expression implies that / tends to zero at x — 0, so that this 
point indeed corresponds to the horizon of the metric. Note the similarity of the expressions 
(51), (53) to the case of the AdS-Schwarzschild metric, Eqs. (26), (24). A careful analysis 



9 As a side remark we note that the system (48) is valid in the case of k ~ H as well. However, the 
boundary conditions (50) take a more complicated form: in this case they must be derived directly from 
the junction conditions (10) and the difference of the position of the brane from the AdS boundary must be 
taken into account. 
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of the obtained bulk metric shows that it indeed describes an AdSs black hole with growing 
mass, as anticipated by the qualitative arguments of Sec. 3.3. 

Once the bulk metric is found one has to substitute it into the modified Friedman equation 
(13): this determines self-consistently the Hubble parameter of the brane. From Eqs. (51), 
(53) we find the next-to- leading order asymptotics of the functions (3 and / at x — *■ 1, 

2 A . flL ^4 



13= 1 + ^(1 -x 

H 1 - x V 80 v 



[l-x) 2 \ 80 

Substitution of these expressions into Eqs. (47) yields that the functions r and F have the 
form (12) in the vicinity of the brane with 



5 \tj 10 \t 

Then from Eq. (13) we obtain the modified Friedman equation (40). 

The following comment is in order. So far in this section we have assumed that the brane 
is completely empty, so that its expansion is purely inflationary. However, the above analysis 
is also applicable if there is matter on the brane provided that the matter density is small, 
G N p < H 2 . 



5 AdS/CFT interpretation 

To understand better the mechanism of self-acceleration it is useful to consider it from the 
point of view of the AdS/CFT interpretation of the RS model [24]. In this framework the 
5-dimensional setup considered above is dual to the 4-dimensional theory described by the 
action (cf. Eq. (4)) 

Sm = J <^x^-g (" I g^ + C-sm + £hs[A;, ...] + C CFT [A^ . . ]\ , 

where Cqft represents the Lagrangian of a strongly interacting conformal field theory cou- 
pled to the vector fields. The properties of this CFT are determined by the structure of the 
bulk sector of the 5d theory. The expression for the energy density of the CFT plasma with 
temperature T has the following form 10 [25], 

37T 3 

PCFT = 8GJ^ T4 ■ (55) 



10 The expression (55) differs by a factor 2 from that of Ref. [25]. This is related to the fact that in contrast 
to [25] where the brane is considered as the "cnd-of-thc-universe" , we assume the brane to be embedded 
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The presence of the [£7(1)] 3 gauge symmetry in the five-dimensional Lagrangian implies that 
the CFT is coupled via [£/(l)] 3 gauge interactions to the vector fields A a . The latter are 
identified with the values of the vector fields on the brane in the 5d setup. 

In the 4d language the condition (5) leads to the following cosmological Ansatz for the 
vector fields: Aq — 0, Af — Ma(t)Sf. This configuration corresponds to three orthogonal 
homogeneous "electric" fields 11 ; in the local inertial frame they have the form, 

E° = MH51 . (56) 

These electric fields give rise to electric currents j a in the CFT plasma. In the linear response 
approximation the currents read 

■f = Cl E a + c 2 TE a , (57) 

where c\, c 2 are dimensionless coefficients. The expression (57) is the most general one 
allowed by dimensional considerations. The coefficients c\, c 2 can be extracted from the 
study of the vector fields in AdS, see appendix C. The result is 

f -glutei , ^>T 

ci = { 7 ; ' (58a) 

C2 = ^k ' (58b) 

where Q is the characteristic frequency of the electric field. For the case of cosmology we 
have Q ~ H <C T and the first term on r.h.s. of Eq. (57) can be omitted. 

The currents heat the plasma with the energy release per unit time per unit volume given 
by 

W = J2 J a E a • ( 59 ) 

a 

Thus, the energy balance equation for the CFT plasma has the form 

Pcft + 4:Hp C FT = W . (60) 

Substituting expressions (55), (59), (56) into Eq. (60), and solving the resulting differential 
equation for the temperature of the plasma we obtain 

1 fSXk 2 f , 2 X 1/3 

T = - I — I dr aa 

Z7T d 



into an infinite space-time consisting of two copies of AdS related by the Z^ symmetry. This leads to the 
doubling of the CFT degrees of freedom as compared to [25] . 

11 To avoid confusion we stress that these are the fields of the hidden sector and have nothing to do with 
the electromagnetic field of the Standard Model. 
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where A is defined in (38). One observes that the contribution 

_ 8nG N 

— ~ PCFT 

of the CFT plasma into the 4-dimensional Friedman equation coincides with the expres- 
sion (39b). 

The analysis of this section suggests a transparent interpretation of the self-acceleration 
mechanism in our model. It is based on the inflow of energy into the conformal matter which 
compensates for cooling of the plasma due to the cosmological expansion. Ultimately, this 
leads to stabilization of the plasma temperature. Similar ideas were expressed previously 
in Ref . [20] ; our model provides explicit realization of these ideas in the framework of field 
theory. Note that the CFT picture also indicates that Eqs. (39) are more robust than could 
be inferred from the 5-dimensional analysis of Sec. 3. Indeed, as discussed in Sec. 3, the 5d 
analysis is not completely reliable as it does not self-consistently take into account the time 
evolution of the bulk metric. On the other hand, the formulas of the present section are 
expected to be valid once the CFT plasma is in thermal equilibrium. This requirement is 
satisfied if the Hubble time 1/ H is much larger than the thermalization time of the plasma. 
By dimensionality, the latter is of order 1/T. As 1/T <C 1/H, we see that the condition of 
thermal equilibrium is always satisfied. 

In our model the inflow of energy into the plasma is produced by the work done by 
electric fields. The key point here is to find a way to prevent the electric fields from rapid 
decay which is usually caused by the cosmological expansion. The condition (5) is crucial 
for this purpose. 

This consideration reveals the following property: the hidden sector responsible for the 
condition (5) has to violate null energy condition (NEC). Indeed, the equation of the energy 
balance for the hidden sector implies 

Phs + 3H(p HS + Phs) = ~W , 

where W is given by Eq. (59). As W > and p^s ~ y at r — > +oo (cf. condition (2) of 
Sec. 2) the NEC is violated: 

Phs+Phs <0 . (61) 

Generically, violation of NEC leads to instabilities [21], so the condition (61) is rather wor- 
risome. We discuss this issue in more detail in the next section. 

Let us point out possible phenomenological generalizations of the self-acceleration mech- 
anism proposed in this paper. One can consider a class of cosmic fluids interacting with 
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electric fields (56) with a general dependence of the conductivity on the density of the fluid, 

j a = c( Pf )E a . 

The case considered in this paper corresponds to c oc pi . The equation of the energy 
balance of the fluid takes the form 

p f + 3H (p f + Pf ) = 3M 2 H 2 c(p f ) . (62) 

Again, one assumes that the constant electric fields (56) are produced by some NEC- violating 
hidden sector, whose energy density pns is zero ( or tends to zero during the cosmological 
expansion). Such a setup will lead to self-accelerated cosmology if the system consisting of 
Eq. (62) and the Friedman equation has an attractor fixed point with pj ^ 0. We point out, 
however, that we are not aware of a microscopic theory of a fluid with c(pf) different from 

1/4 
C OC Pf . 

6 Discussion 

In this paper we proposed a model realizing the self-accelerated cosmology. The model is 
formulated in the context of a brane-world scenario with warped extra dimensions. We 
have shown that the regime of accelerated cosmological expansion is an attractor in the case 
when there is a sufficient energy exchange between the expanding brane and the bulk. This 
exchange is realized by means of bulk vector fields which have a specific interaction with 
the brane. The key requirement for the existence of the self-accelerated solution is that this 
interaction produces the boundary condition (5), which relates the vector fields to the scale 
factor. The condition (5) is analogous to the conditions appearing in effective a-models of 
particle physics. 

The condition (5) is far from being trivial: it breaks explicitly the (Abelian) gauge 
symmetries of the bulk action; moreover, it leads to spontaneous breaking of the Lorentz 
symmetry. Per se, these properties do not imply any inconsistencies. In particular, there 
exist a number of self-consistent models [26, 27, 28] involving spontaneous Lorentz symmetry 
breaking. Still, the stability requirements impose tight constraints on the structure on these 
models. In our case construction of a suitable brane Lagrangian Chs is complicated by the 
fact that it has to produce an energy-momentum tensor violating NEC. 

Let us demonstrate that, under quite generic assumptions, the latter requirement is 
unavoidable in brane-world models with self-accelerated cosmology. Indeed, consider a brane 
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which performs de Sitter expansion in a higher- dimensional bulk. The energy conservation 
equation for the brane has the form, 

p b + 3H(p b + Pb ) = -W , 

where p b , p b are the density and pressure of the matter on the brane and W is the energy 
flux radiated by the brane into the bulk. Let us assume that the brane energy-momentum 
tensor respects the symmetries of de Sitter space-time. This implies p\ = 0. Another natural 
assumption is that W > 0: the brane radiates rather than absorbs energy from the bulk. 
Thus we obtain 

Pb+Pb<0- (63) 

The equality can be attained only if W = 0, i.e. the brane-bulk energy exchange is absent. 
A self-accelerated cosmology of this type is not excluded: an example is provided by the 
DGP model [14]. However, generically one expects W > (this is the case in our model). 
Then, the inequality in (63) is strong and the brane energy-momentum tensor violates NEC. 

It was shown in Ref. [21] that under broad assumptions violation of NEC leads to expo- 
nential instabilities. At first sight, this statement precludes any hope to construct a viable 
brane Lagrangian fitting into the self-acceleration mechanism proposed in this paper. How- 
ever, important assumptions of [21] are that the NEC violating energy-momentum tensor 
is associated with dynamical fields and that one can neglect higher derivative terms in the 
action of these fields. This suggests two possible ways to avoid instabilities. First, one can 
envisage that higher derivative terms in the action for unstable modes cut off the insta- 
bilities at short wavelengths [29] and make them slow enough to satisfy phenomenological 
constraints. Second, some or all of the fields localized on the brane may be constrained 
to vanish by additional requirements which should be attributed to UV-completion of the 
theory. In appendix A we present a model brane Lagrangian exploiting the latter possibil- 
ity. For this Lagrangian we have been able to establish, in a number of phenomenologically 
relevant regimes, the absence of instabilities predicted by [21]. The details of this stability 
analysis will be published in a separate article [22]. 

The model considered in this paper requires UV-completion above the energy scale M ~ 
1 TeV. The latter scale turns out to be the same as the natural scale for the new physics 
stabilizing the Standard Model Higgs mass; so, it is tempting to speculate that the two scales 
could have common physical origin. In this case our seesaw formula (3) for the effective 
energy density corresponding to the late-time cosmological expansion would look appealing 
in light of the cosmic coincidence problem - - the approximate equality of the densities of 
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dark energy and dark matter at the present epoch. Indeed, as shown in [30], the cosmic 
coincidence is a natural consequence of Eq. (3) in the case when the dark matter consists of 
massive particles whose masses and interactions are set by the scale M. Particles with the 
latter properties are generic in the UV extensions of the Standard Model. 

It would be interesting to embed our setup in some higher dimensional (super)gravity 
theory. The following considerations suggest that this may be possible. The [?7(1)] 3 gauge 
symmetry of the bulk action may be interpreted as the isometry of the internal space in the 
compactification of 8-dimensional gravity on a 3-dimensional torus. The bulk vector fields 
then naturally appear from the components of the 8-dimensional metric. The breaking of the 
U(l) symmetries on the brane is also natural and corresponds to breaking of the translational 
symmetry of the internal torus by the position of the brane in it. In this context the condition 
(5) should appear as the consequence of nontrivial winding of the 3-brane around the internal 
torus. We leave the development of this idea for future studies. 

We have considered interpretation of our model in terms of AdS/CFT correspondence. 
The CFT interpretation of the self-acceleration mechanism proposed in this paper allows to 
formulate its phenomenological generalizations. The latter involve consideration of a class of 
cosmic fluids interacting with "electric" fields. It would be interesting to work out predictions 
of this class of models and confront them with observations. 

As a by-product of our study we found an explicit solution of the Einstein equations 
describing a black hole with time-dependent mass in AdSs. We used this solution to embed 
a brane into it. However, this solution can be considered on its own without introduction of 
the brane. In this case the solution is dual by AdS/CFT correspondence to a time-dependent 
configuration of CFT in flat 4d space-time. One finds that it provides a self-consistent 
description of the response of CFT plasma to a triplet of external time- dependent electric 
fields 12 

E a {t) oc 1/t 2 . 

It would be interesting to understand whether this solution can be useful for the description 
of collective phenomena in strong coupling regime of gauge theories. We leave this question 
for future. 
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A An example of the brane sector Lagrangian 

As a possible choice for the Lagrangian of the brane sector of our model we propose 

Chs = XiiQ^d^d^ - 1) + x a \r v AlAl + M 2 5 ab ) , (64) 

where 

and Xi, x ab are Lagrange multipliers. Note that the Lagrange multiplier xi enforces the 
gradient of the field to be a unit time-like vector. This implies that the field <fi is non- 
dynamical; it plays the role of "cosmic time" . 

Let us demonstrate that the hidden sector (64) satisfies the conditions (1), (2) of Sec. 2. 
In the cosmological context one makes an Ansatz 

= 1, d i( p = o. 

Combining with Eqs. (8), (9) we see that in the cosmological setting A a = A a . Then, 
the Lagrange multiplier x ab enforces the condition (1). From (64) one obtains the energy- 
momentum tensor of the hidden sector, 

T^ = 1 Xl d^dA + 2 x ab Affi . 

For the cosmological Ansatz 

Phs = 2%i . 

The evolution of xi is determined from the equation obtained by varying (64) with respect 

tO 0, 

W(xid^-X ab A^A b »d u <j ) )=0. 
For the cosmological setup this equation simplifies and yields 

1 

Xi oc — . 
a 6 
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Thus the condition (2) is also satisfied. 

As discussed in the main text the energy-momentum tensor of the hidden sector violates 
NEC on the self-accelerated solution, so one has to worry about possible instability of this 
background due to rapidly growing short wavelength modes [21]. We have carried out an 
analysis of stability of the self-accelerated background in the model with the hidden sector 
Lagrangian (64); this analysis will be presented in Ref. [22]. The preliminary results of this 
analysis are that the instabilities are absent in the range of wavelengths/time-scales 

(ffM„r 1/2 <A/, Ar^H- 1 . (65) 

We also found that in this range of distances/time-intervals the linearized gravitational field 
of external sources is described by the standard formulas of the four-dimensional Einstein's 
gravity. 

The range (65) of distances/time-scales, where we were able to establish the absence of 
instabilities, is limited both on short- and long-scale sides. This is due to technical difficulties 
of the analysis and does not imply that instabilities appear outside this range. On the other 
hand, for the moment, we are not able to completely eliminate this possibility. 

B Self-accelerated solution: technicalities 

Here we construct the solution of the system of ordinary differential equations (48) with 
boundary conditions (49), (50). First, we consider the system (48) in the interval x G [1— e, 1], 
where e is some small number, e < 1, to be specified later. In this region one can remove 
the explicit dependence on x in the equations by setting x — 1. Multiplying Eq. (48a) by a 
factor of 2 and subtracting it from Eq. (48c) one obtains 

/'Y /7?' + 2^ _ 4 (^ Y /«; Y aa' f : 2[3> 



fj fP P \PJ \PJ P 2 f P 

Divided by /3, this equation can be integrated as 

/' P' „ [ , ( K) 2 ««' /' 2/5' 



where C is the integration constant to be determined from the boundary conditions. We 
will shortly see that the terms in the integrand in Eq. (66) are of the order one or smaller. 
As x changes in the narrow interval of the length e< 1 the whole integral is small and can 
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be neglected. Then, using the asymptotics (50) we obtain, 

J_ + 2^— = 2. (67) 

f/3 (3 2 { ' 

Consider now Eq. (48e). It can be integrated with the result 

( a' + -) = ^ - [ dx 2(3a' = ^ + 0((3a'e) , (68) 

V 2/ a J a 

where fa is an integration constant and we explicitly wrote down the corrections to keep 
them under control. Substituting the expressions (67), (68) into the difference of equations 
(48b) and (48a) we obtain 

where in the last expression we singled out the value of the function a on the boundary and 
accounted for the corresponding corrections. Neglecting the corrections in Eq. (69) we can 
integrate it with the result (51). Then, from Eq. (48d) we obtain the formula (53) up to 
relative corrections of order 0(1/(3). Now, one can check that the equations (67), (48a) are 
satisfied, the latter up to corrections 0((a'/(3) 2 ). Let us consider the asymptotic behavior 
of the metric functions when j3 ~ [3h- From Eqs. (51), (53) we obtain, 

(3 - fa oc fte-^ 1 "*) , / oc /^ e -«-*) . (70) 

We see that the functions (3 and / exponentially approach their asymptotics at 1 —x ^> l/fa. 
We will shortly see that fa ^ 1- Consequently, the parameter e introduced above can be 
chosen to satisfy 1/ fa CeC 1. Finally, from Eq. (68) we obtain at 1/ fa < 1-iCf 

a ' « !% - ££ , a « a , (71) 

where we neglected corrections of order 0(a'e). 

We now solve the system (48) for the values of x which are well separated from the 
boundary, 1 — x ^> 1/ fa. In this region / <^ 1. One sets / = in Eq. (48d) and integrates 
with the result 

/ r \ 1/3 

One chooses the solution which is regular at x — 0. Thus, we obtain that (3 is simply 
constant, (3 = fa. Inserting this result into Eq. (48e) we obtain 
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Again, we should take the constant solution, a = ao- Substituting a' — into Eq. (71) we 
obtain the relation (52). Thus, fa is indeed large, fa >> 1. Equation (48a) is now trivially 
satisfied while Eqs. (48b), (48c) reduce to 



/' a 



h 



f m 

After integration this gives Eq. (54); the latter matches with Eq. (70) provided the relation 
(52) is satisfied. 

Now, one can go over the derivation and check that the corrections which we neglected 
in various expressions are small. For example, consider the corrections in Eq. (69). From 
the relation (52) it follows that a' < a ; taking into account that (5 > fa one obtains 
aa'e/p 2 < (3. 

We have checked the analysis described in this appendix by numerical solution of the 
system (48) with boundary conditions (49), (50). 

C Conductivity of CFT matter 

In this appendix we calculate the conductivity of the CFT plasma which is dual to the 
bulk sector of our model. A similar calculation for the M = 4 supersymmetric Yang-Mills 
plasma was performed in Ref. [31]. First, let us relate the temperature of the plasma to the 
properties of the bulk geometry. By the standard rules of the AdS/CFT duality (see, e.g., 
Ref. [23]) the finite temperature on the CFT side corresponds to the presence of a black hole 
in the bulk. The temperature of the plasma is then equal to the Hawking temperature of 
the black hole [25], 

T = — 

7rr 

Second, we use the relation (59) to derive the formula for the CFT current. Equating W 
to the energy injected by the brane into the bulk, W = 2T t (where a factor of 2 is due to 
^-symmetry) and using Eq. (16b) we obtain, 



2 



Here we canceled the factor 3 appearing both in Eqs. (16b) and (59) due to the presence of 
3 vector fields. Using the relation 

E = Md t A/a 2 
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between the physical electric field and the vector potential we obtain 

This formula expresses the CFT current in terms of the quantities describing the response 
of the bulk system to the boundary perturbations. 

Let us consider periodic electric field, E = Eoe~ tnr . This corresponds to periodic per- 
turbation of the vector fields on the AdS boundary. Note that the physical frequency fl is 
related to the frequency u in conformal time t by Q = uo/a. The response of the bulk to 
such perturbation was calculated in Sec. 3. Let us first consider the case of small frequency, 
1] <T. It is easy to see that this condition is equivalent to Eq. (29), thus we can use the 
formula (35). Substitution of the latter formula into Eq. (73) yields 



3 = <~i§k ln 
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Taking T <^i k and using r 2 /F f» 1 on the brane one obtains the expressions (58a) (lower 
case) and (58b). In the opposite limit Q ^> T the response of the vector field to the boundary 
perturbation is not affected by the presence of the black hole. In this case we can use Eq. (19) 
for the propagation of the vector field in the AdS metric. Performing the same steps as above 
we obtain 

This reduces to Eq. (58a) (upper case) in the regime Q <^ k. 
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